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The  path-integral  successive-collision  method  is  extended  and  applied  to 
the  set  of  Unear  Boltzmann  equations  for  a  very  general  initial-  and  bound¬ 
ary-value  multi-swarm  problem  allowing  for  volume  production,  reactive 
coUisions,  and  various  partkle-surface  interaction  processes.  The  complete 
formal  solution  is  given,  and  its  convergence  is  inferred  from  a  simple 
physical  argument.  - _ _ _ _ 
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The  linearized  Boltzmann  equation  has  been  widely  used  to  describe  the 
transport  of  particle  swarms  in  prescribed  fields  and  background  media,  such 
as  neutrons  in  fissible  material,1  electrons  in  semiconductors,2  or  electrons  and 
ions  in  weakly  ionized  gases.2'*  These  studies  have  two  basic  goals,  namely  (i) 
calculating  transport  coefficients  (i.e.,  moments  of  the  distribution  functions), 
and/or  (ii)  finding  the  distribution  functions  themselves. 

The  present  work  is  concerned  with  the  second  of  these  goals.  We  first 
formulate,  in  integro-dilferential  form,  a  very  general  bounded  and  dynamic 
multi-swarm  problem  which,  to  our  knowledge,  includes  all  situations  consid¬ 
ered  previously  as  special  cases.  The  relevant  linear  Boltzmann  equations  are 
then  cast  into  their  path-integral  forms,  and  their  complete  formal  solution 
is  derived  by  means  of  an  appropriately  extended  version  of  the  successive- 
collision  iteration  method.2,4'2  This  scheme  is  readily  interpreted  in  physical 
terms,  and  its  convergence  is  inferred  from  a  simple  physical  argument. 

Consider  a  fixed  spatial  domain  B  surrounded  by  the  surface  A  as  shown 
in  Fig.  1,  and  assume  that  confined  therein  are  n*  particle  swarms  numbered 
by  the  swarm  (or  species)  index  9  (9  ■»  I, . . . ,  n»).  Let  species  <r  be  character¬ 
ized  by  the  particle  mass  m#,  the  electric  particle  charge  q",  and  the  velocity 
distribution  function  /#(r,v,t).  The  swarm  particles  move  in  prescribed  elec¬ 
tromagnetic  and  gravitational  fields  E(r,  t),  B(r,  t),  and  g(r,t),  at  the  same  time 
undergoing  binary  collisioas  with  the  molecules  of  a  background  gas  having  an 
arbitrary  but  given  velocity  distribution  function.  Allowing  for  volume  produc¬ 
tion  of  swarm  particles,  the  integro-dilferential  linear  Boltzmann  equation  may 
be  written 
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where 

*'(r,v,t)  -  -£•  [e(M)  +  -  X  B(r,t)l  +  g(r.t) 

(B  L  C  J 

is  the  acceleration,  the  source  term  7*(r,v,t)  describes  volume  production, 

G"'(r,v,t)-  j  <?Y'K"'(r,Y+-v',t)r'(ry,t) 

describes  the  gain  for  species  *  from  collisions  of  species  o' ,  and  P*(r,v,  t)  is 
the  total  collision  frequency.  The  scattering  kernel  v  ♦—  v*,  t)  basically 

represents  the  probability  for  a  species-*1  particle  with  pre-collision  velocity  V 
to  end  up  after  collision  as  a  species-*  particle  with  velocity  v.  Any  kind  of 
collision,  conservative  or  reactive,  is  admitted  for  which  the  GTs  and  P’s  remain 
Unite. 

Equation  (1)  must  be  complemented  by  initial  and  boundary  conditions. 
Let  us  assume  that  at  the  initial  time  the  swarm  distributions  are  given 
everywhere  as 

f*[r,Ytt «  tf) »  /f  (r,  v) ,  (2) 


and  that  the  relation  between  incoming  and  outgoing  particles  at  the  surface  A 
is  of  the  fairly  general  form 
/a(p  a,  v<«,  t)  —  g\(rA,  Y(n,  t) 

+  £  f  {r  a, Yin*- Y,t%i,t)f*'{rA, 


where  r*  is  a  “surface  point”  (i.e.,  a  position  still  within  B  but  infinitesimally 
close  to  A),  fA{rA,Yin,t)  a  fr{rA,Y(%,t)  (which  notation  is  convenient  because 
!A{rA,Ytn,t)  will  be  treated  as  an  unknown),  and  v<»  and  v0><  are  inward 
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and  outward  directed  velocities,  respectively,  cf.  Fig.  1.  The  (given)  func¬ 
tion  jJIrx.Vf*,*)  describes  externally  controlled  injection  into  B,  while  the 
“generalized-  reflection”  kernel  (r a.  v<»  «-  v£,t,i)  basically  represents  the 

probability  for  a  species-^  particle  hitting  the  surface  A  with  velocity  to 
be  “re-injected”  as  a  species-9  particle  with  velocity  v<».  Clearly,  this  model  is 
capable  of  describing  such  processes  as  elastic  or  inelastic  reflection,  secondary 
emission,  surface  ionization  and  neutralization,  or  simply  absorption  of  swarm 
particles. 


Equations  (1),  (2)  and  (3)  define  a  very  general  evolution  problem  for  the 
functions  ▼<*,*)  and  /'(r,v,t),  which  we  now  cast  into  a  path-integral 

form  suitable  for  both  general  considerations  and  practical  applications.  .Vs  is 
well  known,  the  operator 


dt 


8 

—  +  v  ■ 
dt 


(4) 


appearing  in  (1)  represents  the  time  derivative  associated  with  the  motion  of 
species-^  particles  along  their  collisionless  phase  trajectories  (free  paths,  charac¬ 
teristics).  The  species-9  characteristic  through  the  point  [r,v,t|,  to  be  denoted 
by  (r*(t),  v*(t),t],  can,  in  principle,  always  be  calculated  from  the  equations  of 
motion  dt*/dt  ■■  v#  and  iv*  jdi  *■  ir,  with  t  <  t  and  the  “Anal"  conditions 
—  t)  r,  v<r(t  v  e)  »  ▼.  Of  coarse,  the  symbol  a'  introduced  here  is 
shorthand  for  a'(r#,  r#,  t),  and  analogous  notations  will  be  used  for  other  quan¬ 
tities,  too.  With  (4),  formal  integration  of  (1)  along  the  species-9  characteristic 
through  |r,  v,  tj  leads  to 


/*(r,v,«)  —  /*(*?,  exp(— *'(fg,<)| 


exp(-«#(<,t)| , 


(5) 


f 


-5- 


where  [r*  —  r*(t*),  **  "  —  t*]  is  the  (suitably  chosen)  “origin"  of 

the  species**  characteristic  through  (r,  v,  t],  and 

*'(*!,*)-  f*  dt9* [r*,v*,i) 

Jit 

is  the  “collisional  relaxation  exponent”.  In  some  eases  —  e.g.,  in  time-independ¬ 
ent  problems — it  may  be  convenient  to  re-write  the  path  integrals  by  using  one 
of  the  identities  it »  |dfr*|  /  |V'|  and  dt  —  |dv#|  / 14'|. 

Equation  (5)  represents  the  path-integral  version  of  (1)  and  is  readily  inter¬ 
preted  in  physical  terms.  As  one  proceeds  along  a  characteristic,  the  memory 
of  conditions  at  the  origin  is  gradually  wiped  out  by  collisional  loss  (first  term 
on  the  r.h.s.),  while  a  new  distribution  builds  up  through  both  collisional  gain 
and  volume  production,  with  more  recent  contributions  being  dominant  (second 
term). 

The  species-*  characteristic  passing  through  the  point  [r,  v,  t]  can  either 
originate  at  an  “initial”  point  [r*  a  r*(t  *■  t,),  vf  a  v'(t  —  t,),t  —  t4],  with  rf 
lying  in  the  interior  of  B,  or  at  a  surface  point  [rj  a  r#(t  ■»  t\),  rf,  a  v#(t  ■■ 
t^),i  *  t\\,  with  tf  <  <  t.  These  trnjectories,  and  the  points  through 

which  they  pass,  will  be  referred  to  as  “i-type"  and  “A-type",  respectively,  cf. 
the  curves  C\  and  C*  in  Fig.  1.  For  what  follows  it  is  convenient  to  introduce 
the  unit-step  functions  C/*(r,r,  t)  (equals  unity  for  i-type  points,  vanishes  for 
A-type  points)  and  £/JJ(r,  v,  t )  (equals  unity  for  A-type  points,  vanishes  for  i-type 
points). 

With  these  preparations,  Eq.  ($)  can  be  re-written  in  more  detail  as  follows: 


-  - 
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r  -  Uf  |/T  exp[— t)|  +  f'd  G”'  +  7'  exp|-*'(M)|  J 

+  UA  | /a  expl-ie'lt^.t)!  +  dt  ^  +  7*  «*P|-*<r(M)|  j 


(#) 


where  /*  —  /®(rf,  vf,t<)  and  j\  —  /#(*,5«^i5 .O-  This  integral  equation 
is  very  convenient  in  that  it  explicitly  accounts  for  the  mixed  initial-  and 
boundary-value  character  of  the  problem.  It  has  to  be  solved  simultaneously 
with  the  boundary  condition  (3),  which  also  becomes  an  integral  equation  for 
/5(*>t.v<«.f)  and  /#(r,  v,t)  when  /*'(r/t  « vo*t>0  °n  the  r.h.s.  is  re-expressed 
in  terms  of  these  functions  via  (8). 


We  now  solve  Eqs.  (3)  and  (8)  formally  using  the  method  of  successive 
collisions.2-4"7  To  this  end  we  decompose  the  two  unknown  distribution  func¬ 
tions  in  the  form 


00 


Ia{*a, »<*,«)  -  £ 

t 

(7a) 

00 

f*(r,Y,t)  -  £/?(r,r,f) 

(76) 

where  the  "partial  distribution  functions"  f\  v  and  /'  are  the  distribution  func¬ 
tions  of  the  “cla*s-(*,  i/)"  particles,  by  which  we  mean  those  species-*  particles 
that  have  suffered  exactly  *  generalised  collisions  since  their  first  appearance 
in  the  system.  By  definition,  a  swarm  particle  suffers  a  “generalized  collision' 
either  when  it  collides  with  a  background-gas  molecule,  or  when  it  hits  the  sur¬ 
face  A.  By  “first  appearance"  of  a  swarm  particle  we  mean  either  its  presence 
within  f  at  t  «  If,  its  creation  by  volume  production  at  some  I  >  tt,  or  its 
injection  through  the  surface  at  some  t  >  t«.  Classes  with  v  «  0  are  populated 
by  these  first-appearance  mechanisms,  whereas  classes  with  *  >  1  gain  their 


particles  from  classes  with  v  —  I  through  generalized  collision  processes.  Ev¬ 
ery  class  (<r,  v)  loses  particles  through  generalized  collision  processes,  and  if  the 
latter  are  non-destructive,  these  particles  go  into  classes  with  v  +  i.  According 
to  this  philosophy,  the  partial  distribution  functions  must  clearly  satisfy  the 


equations 


/A.olrA,  t )  —  gJ(PA.  Vi,,  f  ) 

/»  —  Vi  {/<  «P(-*'(*(.*)|  +  ^  *7'  exp(— 

+  Ua  {&  exp(-«'(£,t)|  +  <*7*  exp(-*'(M)l  j 


for  v«0,  and 


v„  —  v;.t, t)/zh (p^  >ve*r>*)  (®a) 
ft  *  [  X!  ^v-i]  «xp(-*'(M)| 

*'  Lr'"‘  •*  (96) 

+  &a  j Ta.u  +  Jir  &  Ol-i  exPl-K'MjJ 

for  u  >  1,  where 

Gl'h  -  /  dV  fC"'(r,  v  -  v',  t)  /'Ijr,  v',  t) .  (10) 

Equations  (8a,b)  yield  explicit  expressions  for  Saa  /•  >  whence  G%9'  can 
be  calculated  via  (10).  Using  these  results  in  (9a, b)  and  (10),  one  successively 
finds  /JJ  /*  and  Gi0>,  which  in  turn  can  be  inserted  into  (9a,b)  and  (10) 
to  yield  /a. j,  /j  and  Gf#\  etc.  Thus,  Eqs.  (7)— (10)  define  a  closed  iteration 
scheme,  which  yields  the  formal  solution  to  our  swarm  problem  and  can  be  seen 
to  converge  as  follows. 
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In  a  specific  application,  the  number  of  partial  distribution  functions  that 
must  be  retained  in  (7a.b)  for  a  meaningful  approximation  to  the  exact  distri¬ 
bution  functions  increases  with  the  number  of  generalized  collisions  a  typical 
swarm  particle  can  undergo.  E.g.,  in  a  situation  where  there  are  hardly  any 
swarm  particles  having  undergone  more  than  three  generalized  collisions,  one 
may  safely  neglect  in  (7a,b)  ail  /£  „’ s  and  /*’s  with  say  v  >  i.  This  indicates 
that  the  method  of  successive  collisions  may  become  computationally  expen¬ 
sive  in  problems  where  the  particles  have  to  be  traced  through  many  collisions. 
However,  for  finite  evolution  times  (t  -  U  <  oo)  the  scheme  (7)-{10)  is  always 
bound  to  converge,  and  since  it  can  handle  non-standard  (e.g.,  strongly  non¬ 
equilibrium)  situations  in  a  natural  way,'-7  it  may  well  provide  a  suitable  basis 
for  attacking  swarm  problems  which  because  of  their  complexity  have  hardly 
been  looked  at  as  yet.  As  an  example  we  may  mention  swarm  behavior  near 
walls,  which  is  important  but  needs  further  consideration.* 

In  conclusion,  recent  studies*-7  and  the  present  work  have  contributed  to¬ 
wards  making  the  path-integral  successive-collision  method  a  practicable  tool 
of  transport  theory,  which  seems  suitable  for  attacking  more  complex  swarm 
problems  than  hitherto  considered. 
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